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Abstract 

We compute the leading four physical terms in the low-energy expansions of 
heavy-light quark current correlators at four-loop order. As a by-product we re¬ 
produce the corresponding top-induced non-singlet correction to the electroweak 
p parameter. 
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1. Introduction 

Two-point correlation functions of heavy-light quark currents have found 
use in a number of phenomenological applications. One example is the predic¬ 
tion of corrections to the electroweak p parameter a a a , where the flavour 
non-diagonal correlator of vector currents is required for vanishing external mo¬ 
mentum. Another important class of applications is the sum-rule determination 
of meson decay constants (see e.g. [4|, |5|). Here, the absorptive part of the 
respective correlators above the production threshold is needed. 

Progress in lattice simulation may allow precise determinations of even more 
QCD parameters. For instance, the values of the strong coupling constant, 
the charm quark mass and the bottom quark mass have been determined with 
high accuracy from moments of heavy-heavy correlators in E0. In these 
analyses, moments of flavour diagonal currents have been determined on the 
lattice choosing a frame where the spatial momentum of the correlators vanishes. 
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The values of the quark masses and the coupling constant are then extracted by 
equating these moments to their counterparts calculated in perturbation theory 
at the four-loop order 


nmmEiQ- 


The methodology is thus similar to traditional quarkonium sum rules |15l . 
Tl . 171. but using lattice moments in place of moments of the experimentally 
measured hadronic R ratio. While for the sum rules only the correlator of vector 
currents can be used, there is no such restriction for the lattice simulation. In 
fact, in Q different Lorentz structures were considered, with the most precise 
results stemming from pseudoscalar currents. Furthermore, also correlators of 
heavy-light currents could be used to extract the values of the charm and bottom 
quark masses and possibly the strong coupling constant |l8j. To be competitive 
with the analyses for the heavy-heavy case the corrections to the perturbative 
moments of the heavy-light current correlators have to be known up to four 
loops. These corrections are presented in this work. 

Given their usefulness, perturbative corrections to heavy-light correlators 
have been studied quite intensively and analytic results up to two loops have 
been known for many years [lfj, |20|. While the three-loop correction is not 
known analytically, many terms in expansions in both the low-energy and the 
high-energy limit have been calculated in 21, 22, 23|. Combining these with 
the behaviour near threshold, accurate approximations for arbitrary kinematics 


have been constructed [21], [22|. In the low-energy region also corrections due to 
a non-vanishing light quark mass are known [24l . [251. 

The four-loop corrections remain mostly unknown. In the high-energy region 
the leading term is equal to the non-singlet part of the corresponding diagonal 
correlator, which has been computed for both scalar and vector currents 


l28l |. In the low-energy region, conversely, there is no such simple correspondence 
between diagonal and non-diagonal correlators. The vector correlator in the 
limit of vanishing external momentum constitutes a central ingredient in the 
determination of non-singlet four-loop corrections to the p parameter, which 
have been calculated in [2], y] ■ 

In this work we present the four-loop corrections to the low-energy expan¬ 
sions of both scalar and vector heavy-light quark current correlators up to the 
eighth power of the external momentum. After introducing our conventions in 
section [2] we briefly describe the calculational setup and present our results in 
section[3] Section[4]describes the re-calculation of the top-induced contributions 
to the electroweak p parameter, which constitutes an important consistency 
check. We conclude in section [5] 


2. Conventions 

The correlators of heavy-light vector and scalar currents are defined as 


n M „(g) = i J dx e iqx (0\j^x)j„{Q)\0) 

n {q) = i [ dxe lqx (0\j(x)j(0)\0) 


(1) 

( 2 ) 


2 






















with the vector current j t i(x) = V>(a;)7/iX(0) and the scalar current j(x) = 
We consider a heavy quark if> with the pole mass m and a massless 
light quark y. It should be noted that in the limit of a vanishing light-quark 
mass the correlators of two axial-vector or pseudo-scalar currents coincide with 
the vector and scalar correlators, respectively. 

It is convenient to introduce polarisation functions 

IW q) = (- q 2 g^ + q»q v W(q 2 ) + q^li v L (q 2 ), (3) 

n (q)=q 2 W(q 2 ). (4) 

In the following we will not consider the longitudinal polarisation IV^(q 2 ). The 
perturbative expansions of IT 5 (q 2 ) with 5 = v,s up to four loops read 

irV) = n 5 ’(% 2 )+— 1 1 n 5 ^ 3 )( 9 2 )-i-.... 

7r \ 7T J \ 7T J 

( 5 ) 

Being interested in the limit q 2 —> 0, we can expand the coefficients in the above 
series as 


II 4 ’W(g 2 ) 


3 

167T 2 



3 

167T 2 



( 6 ) 


where we have used the abbreviations z = q 2 /m 2 ,z = q 2 /fh 2 with fh denot¬ 
ing the mass of the heavy quark in the MS scheme. Note that the coefficients 
with n = —1,0 still contain poles in the limit e = (4 — d)/2 —> 0. In physical 
observables these have to be cancelled by the wave-function and mass renormal¬ 
isations of the particles (e.g. W bosons) coupling to the respective current. In 
the following we will describe the calculation of the coefficients Cn ,C n A for 
n < 4. 


3. Calculation and results 


First, the four-loop diagrams contributing to the polarisation functions are 
generated with QGRAF [ 29 }. In the following steps we perform several algebraic 
manipulations with the help of TFORM [30l |31|. As a first simplification, we 
apply partial fractioning to denominators that differ only by their mass and the 
external momentum q, i.e. we use 


1 1 
p 2 (p ± q) 2 — m 2 


__I_(T 

q 2 ± 2 pq — m 2 \p 2 


1 

Jp±qy 



( 7 ) 


Since we will perform an expansion in q the prefactor on the right-hand side 
has no influence on the tadpole topology of the considered diagram. Performing 
partial fractioning before the identification of the diagram topologies greatly 
reduces both the number and the complexity of the topologies that have to be 
considered. Using the algorithm described in Appendix A| we map the resulting 
diagrams onto 28 topologies. 
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Next, colour factors are calculated using the FORM [30] package color [321. 
We choose a routing for the external momentum q which minimizes the number 
of propagators depending on q. After this we evaluate the traces over gamma 
matrices and perform a Taylor expansion in q. The scalar integrals we ob¬ 
tain after tensor reduction and the elimination of reducible scalar products are 
reduced to master integrals using a private implementatiorO] [35| of Laporta’s 
algorithm |36| . All required master integrals are known analytically or numeri¬ 
cally [ 37 I ■ihI. f.iol. 

For the presentation of our results we impose the overall renormalisation con¬ 
dition IT(0) = II s (0) = 0. The corresponding divergent subtraction terms are 
listed in |Appcndix B| For the remaining coefficients according to equation © 
we obtairFI 


c v,(3) = + 14 5508 C 2 Cp + g 4892 CaC 2 + 0 351 C 3 

- 0.22MC A C F T F n h - 0.6242 C 1 2 F T F n h 

- 12.56835 C A C F T F m - 3.07525 C 2 F T F m 

+ 0.107 C F Tjn\ + 0.14 C F Tjn h m + 1.91917 CpTjnl , ( 8 ) 

c v,(3) = + 7.39 iwclCp + 5.65943 C A C 2 F + 0.80504 

+ 0.0683 C A C F T F nh - 0.3114 C 2 F T F n h 

- 6.0806 C A C F T F m - 2.2303 C F T F m 

+ 0.008 C F Tln 2 h - 0.0052 C F Tjn h m + 0.9442 CpTjnf , (9) 

C 3 ,(3) = + 4.42563 C\C F + 3.86666 C A Cj + 0.73105 C% 

+ 0.0448 C A C F T F n h - 0.1713 C 2 F T F n h 

- 3.57037 C A C F T F m - 1.5461 CjTpm 

+ 0.0017 C F Tjn 2 h - 0.005 C F Tjn h m + 0.56396 CpTjnf , (10) 

c v,(3) = + 2.90512 C\C F + 2.7515 C A C 2 F + 0.5965 Cj 

+ 0.0278 C A C F T F n h - 0.104 C 2 F T F n h 

- 2.31867 C A C F T F m - 1.10502 C 2 F T F m 

+ 0.0006 C F Tjn 2 h - 0.003 C F Tjn h m + 0.3708 CpTjnf , (11) 

C-i’ (3) = + 1.6424 C\C F + 1.65318 C A Cj + 1.41042 Cj 

- 1.39916 C A C F T F n h + 0.5551 C 2 F T F n h 

- 3.129699 C A C F T F m + 0.556834 CjT F m 

+ 0.376 CpTjnl + 0.65308 C F Tjn h m + 0.441495 CpTjnf , (12) 


1 The implementation is written in C++ and uses GiNaC fO and fermat f->4( . 

2 The moments are presented in a form that yields five significant digits for QCD with 
n l < 5 light flavours. All results are attached in electronic form as ancillary files to this 
preprint. 
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C* 2 ,(3) = + 5.66925 C\C F + 5.36995 C A Cj + 2.1099 C% 

- 0M76C A C F T F n h + 0.1338 C 2 F T F n h 

- 5.00716 C A C F T F m - 1.60465 C%T F m 

+ 0.037 C F Tj,n 2 h + 0.0314 C F Tjn h m + 0.711301 CpTjnl , (13) 

C|’ (3) = + 4.18695 C\C F + 4.9201 C A Cj + 2.0783 C% 

+ 0.0196 C A C F T F n h - 0.0103 CjT F n h 

- 3.5077 C A C F T F ni - 1.7089 CjT F m 

+ 0.009 C F Tjn 2 h + 0.0006 C F Tjn h m + 0.5215 CpTjnl , (14) 

C s / 3) = + 2.945114 C\C F + 3.81782 C A Cj + 1.6905 Cj 
+ 0.01994 C A C F T F n h - 0.0347 C 2 F T F n h 

- 2.41729 C A C F T F ni - 1.3927 CjTpm 

+ 0.0034 CpTjnl - 0.0016 C F Tjn h m + 0.36976 CpTjnl , (15) 

where we have set the renormalisation scale /i to the on-shell mass m. We follow 
the usual convention for the colour factors with C A = 3 ,C F = 4/3, T/ = 1/2 for 
QCD. The number of light (massless) quark flavours is denoted by rq, whereas 
nh stands for the number of heavy flavours. 

If we choose to express the polarisation functions in terms of the MS mass 
m at the scale /i = fh and a s (fh ), we arrive at 

Cj {3) = _ 1 2994791 C \Cp + 1.20957 C A Cj + 0.537098 Cj 

- 1.75125 C A C F Tpn h + 1.29201 CjT F n h 
+ 0.530618 C A C F Tpm - 0.0193 C 2 F T F m 

- 0.0853 CpTjnl + 0.07322 C F Tjn h ni - 0.0389 CpTjnl , (16) 

C v 2 m = - 1.0623284 C A C F + 1.035507 C A Cj + 0.1608 C% 

- 0.74336 C A C F T F n h + 0.72663 CjT F n h 
+ 0.905515 C A C F T F m - 0.46186 C 2 F T F m 

- 0.0944796 CpTjnl - 0.04078 C F Tjn h m - 0.10011 CpTjnl , (17) 
C"’ (3) = - 0.8577978 C\C F + 1.1607802 C A Cj - 0.1496539 C% 

- 0.4624948 C A C F T F n h + 0.520167 CjT F n h 
+ 0.7959545 C A C F T F m - 0.631599 C 2 F T F m 

- 0.06236 C F Tjn 2 h - 0.027228 C F Tjn h m - 0.08873276 CpTjnl , 

(18) 

= _ 0J17803 C^Cf + 1.286187 C A Cj - 0.40493 C% 

- 0.320038 C A C F T F n h + 0A1065CjT F n h 
+ 0.67538025 C A C F T F m - 0.72124 C 2 F T F m 

- 0.04337 C F Tjn 2 h - 0.018157 C F Tjn h m - 0.076781 CpTjnl , (19) 
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Ci’ (3) = + 1.642401 C\C F - 0.510074 C A C 2 F + 1.41042 C% 

- 1.39916 C A C F T F nh + 1.34177 C 2 F T F n h 

- 3.1297 C A C F T F m + 1.343472 C 2 F T F m 

+ 0.3759 C F Tj,n 2 h + 0.65308 C F Tjn h m + 0.4415 CpTjnf , (20) 

C°' {3) = + 0.38582 C\C F + 0.31725 C A Cj + 0.916757 Cj 

- 0.55487 C A C F T F n h + 0.80004 C 2 F T F n h 

- 0.640838 C A C F T F m + 0.44938 C 2 F T F m 

- 0.02698 C F Tjn 2 h + 0.0092 C F Tjn h m + 0.05861 CpTjnf , (21) 

C f g’ (3) = - 0.039796 C 2 a C F + 0.48671411 C A C 2 F + 0.4018122 C% 

- 0.38621 C A C F T F n h + 0.46308 C 2 F T F n h 

- 0.014653 C A CpT F m + 0.042405 C 2 F T F m 

- 0.0422 C F Tjn 2 h - 0.01723 C F Tjn h m - 0.00067 CpTjnf , (22) 

C s / 3) = - 0.224943 C 2 A C F + 0.565757 C A C 2 F + 0.067439 C% 

- 0.284427 C A CpT F n h + 0.3319473 C 2 F T F n h 
+ 0.2025025 C A C F T F m - 0.162445 C 2 F T F m 

- 0.035056 C F Tjn 2 h - 0.014915 C F Tjn h m - 0.021858 CpTjnf . 

(23) 


4. The p parameter 


To verify the correctness of our calculation we have performed a number 
of cross checks. Obviously, our results are UV-finite. We have also performed 
an expansion up to linear order in the gauge parameter and verified that the 
gauge dependence cancels in the coefficients Cj l ' 3 \cj < ' 3 K The strongest check, 
however, is the comparison to the known four-loop non-singlet corrections to 
the p parameter [2, Q. 

The electroweak p parameter has been introduced in Ref. (40| . Considering 
only QCD corrections it can be written as 

p = 1 + Sp (24) 


with 


5p = 


IIzz(0) n H A l y(0) 


M 2 




(25) 


where II Z z (0) and Tlww (0) denote the self energies of Z and W boson, respec¬ 
tively. 

In order to calculate the contribution from the Z-boson self energy to the p 
parameter we also need the leading moment of the flavour diagonal correlator. 
To this end we introduce IIj iag (< 7 2 ) similar to Eq. JT|) but with the heavy-heavy 
axial current 


3a = ■ 


(26) 
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and the moments 


tt“.(3) 

il diag 


(q 2 ) 


3 

167T 2 




(27) 


In what follows we will only consider the top-induced four-loop correction to p , 
corresponding to pz in the expansion 


Sp = 3 x t ^2 


i—0 


Xt 


V2G F m 2 t 

167T 2 


(28) 


The corresponding corrections to the Z and W self energies then read 


4 3 k o) 

Ml 


= 3 x t 


, _ ^-,0,(3) _ 1^,(3) 

1 d J 7 M T 


l,diag — l,diag 


+ singlet terms, 


and 


n 


(3) 

WW 


( 0 ) 


Ml 


= 3 x t 


i --Ac’ 


~tV,{ 3) 


_ _j_/-<«,(3) 

d Cb ’- 1 


(29) 


(30) 


where the higher-order corrections U^ Z ,H^ W are defined in analogy to equa¬ 
tion ([5]) . Cl'2\ and C 0 jff\ diag denote the moments with n = —1 of the respective 
longitudinal polarisation functions; from an explicit calculation we obtain 


W3) 
-l,- 1 


-<n,(3) 


C 


a,(3) 

L, — l,diag 


= -c 


o,(3) 

— l,diag ‘ 


(31) 


Note that in the non-diagonal case the vector and axial-vector correlators con- 
incide and that the (—l)-th moment of the diagonal vector correlator vanishes. 
The contributions from W- and Z-boson self energies are divergent on their 
own and only their sum is finite. The singlet terms calculated in Ref. [ 2 ] are 
finite on their own and we do not repeat them here. Using the results given 


in Appendix B we obtain in the MS scheme 


P3,non-singlet = ~ G -l’ = 1-60667, 

and after converting to the on-shell scheme 

P3,non-singlet = 101.083, 

in full agreement with the results in the literature 0 , 1 - 


o,(3) 


=»n,(3) 


(32) 


(33) 


5. Conclusion 

We have calculated the four-loop QCD corrections to the low-energy mo¬ 
ments of flavour non-diagonal current correlators up to n = 4. Our results are 


7 












valid for (axial-)vector and (pseudo-)scalar currents in the limit of a vanish¬ 
ing light-quark mass. As a by-product we have confirmed the results for the 
non-singlet correction to the electroweak p parameter first obtained in [ 111 - 
In combination with lattice simulations, our results can be used for the pre¬ 
cision determination of heavy-quark masses. Furthermore, they can serve as 
an ingredient in the approximate reconstruction of the four-loop corrections for 
arbitrary external momenta. For the latter application, however, more input 
from other kinematic regions is still required. 
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Appendix A. Symmetrisation 


The closely related problems of symmetrisation and mapping diagrams to 
topologies are ubiquitous in multiloop calculations. Commonly used algorithms 
employ either the diagrams’ parametric representations 41] or representations 
as graphs. To avoid cumbersome transformations, we choose to work with the 
original algebraic form obtained directly from the Feynman rules. 

A general L-loop scalar diagram I with P propagators has the form 


1 = J[dh]...[dl L ] 


D * 1 ... Dp P 


(A.l) 


with (not necessarily positive) integers a\,...ap. The [dU] are suitable d- 
dimensional integral measures, e.g. as in equation (IB.II) . and the propagators Di 
are functions of the loop momenta l \,... ,lp, any number of external momenta, 
and a mass rrii. Obviously, / is invariant under a change of variables 


M : Zi i—>■ — Mij lj + qi 


(A.2) 


with |det(M)| = 1 and constant vectors g*. 

Consider now a diagram I with propagators D i ,..., Dp and the diagram / 
as defined by eq. (IA.1I) . Let us denote the propagators we obtain by changing 
the loop momenta in / according to eq. (IA.2I) as D[,... ,D' p . We say that 
I and I belong to the same topology iff there is a transformation A4 such 
that {D ' 1: ..., D' P } = {D i,..., Dp}. Likewise, I belongs to a subtopology of 
I iff for some M we have { D\ ,..., D' P } C {Di, ..., Dp}. The problem of 
mapping a diagram to a topology thus reduces to finding out whether a suitable 
transformation A4 exists. 

The basic idea behind our algorithm is to first look for L propagators Di that 
depend on all loop momenta l±,... ,Il- Then we select L appropriate mutually 
different target propagators Dj i and define A4 such that D[ = Dj i . If the sets of 
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the remaining propagators are also equal after applying A4, the two topologies 
are the same. 

To be more concrete, let us now consider a diagram I defined as in equa¬ 
tion (IA.1I) with propagators of the form Di = pf ± to? , where the pi are linear 
combinations of loop momenta and external momenta. The generalisation to 
other forms of the propagators should be straightforward. In practice, we can 
choose the first L propagators to be of the form Di = If ± to? . The algorithm 
then works as follows. 

1. Select a new target topology and choose a representative with propagators 
{£> 1 ,... ,Dp} of the form Di = pf ± to? from it. 

2. Choose a tuple (Di i: ... ,Di L ) (that was not chosen before) of L distinct 
propagators with compatible masses, i.e. = mi,... ,rhi L = rriL • If 
this is not possible go back to step 1. 

3. Consider the next among the 2 L transformations that map the propagators 

(Di,..., D l ) onto (D il ,...,D iL ), i.e. Ij ±p tj j = 1 If no 

transformation is left go back to step 2. 

4. Apply the current transformation to the propagators D i,..., Dp. I then 
belongs to the current target topology if { D [,..., D' p } = {D \,..., Dp}. 
Else go back to step 3. 

As far as identifying the topology of an integral is concerned the algorithm 
terminates as soon as step 4 is completed successfully. For symmetrisation we 
would skip step 1 and always go back from step 4 to step 3 in order to find all 
automorphisms. 

Appendix B. Subtraction terms 

Since the leading coefficients with n = — 1, 0 in equation © still depend on 
the dimensional regulator e = (4 — d)/2, we first have to specify our renormali¬ 
sation prescriptions in d dimensions in order to give meaningful expressions. 

Our d-dimensional integration measure is given by 



(B.l) 


where 7 p ~ 0.5772157 is the Euler-Mascheroni constant. The counterterms in 
the MS scheme are now defined such that they exactly cancel the poles in e. 
For the sake of simplicity, we refrain from defining on-shell renormalisation and 
present the divergent coefficients in terms of the MS quark mass. Writing 



3 —n — *5, (3) 


(B.2) 
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we obtain for fi = m 

= + 1.740 C\C F - 9.555 C A C 2 F + 15.433 C% 

- 7.803 C A C F T F n h + 7.355 C 2 F T F n h 

- 0.228 C A C F T F m - 1.897 C 2 F T F m 

- 0.935 C F Tj.n 2 h + 0.735 C F Tjn h n, + 1.024 CpTjnl , (B.3) 

C—i A = - 1.196 C\C F + 0.592 C a C 2 f - 1.377 C% 

+ 1.130 C A C F T F nf + 0.015 C%T F nf + 0.009 CpTjn 2 , (B.4) 

= + 2.195 C\C F + 0.649 C A C 2 F + 1.278 C% 

- 1.623 C A C F T F n f - 0.244 C 2 F T F n s - 0.025 CpTjnj , (B.5) 

cj’i (3 J = - 1.058 C\C F - 1.750 C a C 2 f - 0.352 C% 

+ 0.635 C A C F T F n f + 0.531 CjTpnf - 0.069 CpTjn) , (B.6) 

= + 0.210 C\C F + 0.516 C a C 2 f + 0.281 C% 

- 0.153 C A C F T F n f - 0.188 C 2 F T F n f + 0.028 CpTjn) , (B.7) 

Co,o 3) = - 0.832 c aPf - 3.606 C a C 2 f + 2.628 C% 

- 1.432 C A C F T F n h + 2.335 C 2 F T F n h 

+ 2.239 C A C F T F m + 0.666 C 2 F T F m 

- 0.425 C F Tjn 2 h - 0.479 C F Tjn h m - 0.330 CpTjnf , (B.8) 

c£’} 3) = + 0.277 C\C F + 0.065 C A C 2 F - 0.180 C% 

- 0.417 C A C F T F nf + 0.172 CjTpnf - 0.020 C F Tjn) , (B.9) 

Co/ 2 (3) = - 0.230 C\C F + 0.019 C A C 2 F 

+ 0.150 C A C F T F n f + 0.024 CjT F n f - 0.017 C F Tjn) , (B.10) 

Co/ 3 (3) = + 0.070 C\C F - 0.051 C A C F T F n f + 0.009 C F Tjn 2 , (B.ll) 

c®’^ = - 72.707 C\C F - 114.585 C A Cj + 20.766 C% 

+ 14.819 C A C F T F n h + 101.776 C 2 F T F n h 

+ 62.816 C A C F T F m + 19.095 C 2 F T F m 

- 17.829 CpTjnl - 24.041 C F Tjn h m - 3.175 CpTjnf , (B.12) 

c!i (3) = - 5.959 C\C F + 10.188 C A C 2 F - 6.959 C% 

+ 16.536 C A C F T F n h - 0.578 C 2 F T F n h 

+ 2.295 C A C F T F m + 0.422 C 2 F T F m 

- 0.544 C F Tj.n 2 h - 0.644 C F Tjn h m - 0.100 CpTjnl , (B.13) 
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cj[ 3 ) 2 = + 7.939 C 2 A C F + 1.310 C A C 2 F + 3.731 C% 

- 7.673 C A C F T F n h - 3.327 C%T F n h 

- 5.840 C A C F T F m - 0.327 C F T F m 

+ 0.481 C F Tln 2 h + 0.296 C F Tjn h m - 0.185 CpTjnl , (B.14) 

= - 3.813 C\C F - 11.708 C A Cj - 2.438 Cj 

+ 2.236 C A C F T F nf + 3.042 C 2 F T F n s - 0.222 CpTjn) , (B.15) 

c*’i ( f2 = + 0.840 C\C F + 4.125 C A C 2 F + 4.500 C% 

- 0.611 C A C F T F nf - 1.500 C 2 F T F n f + 0.111 CpTjn) , (B.16) 

c 0 S o 3) = - 1.740 C\C F + 9.555 C A Cj - 15.433 C% 

+ 7.803 C A CpT F n h - 7.355 C 2 F T F n h 

+ 0.228 C A C F T F m + 1.897 CjT F m 

+ 0.935 C F Tjn 2 h - 0.735 C F Tjn h m - 1.024 CpTjnl , (B.17) 

C 0 S ; 1 (3) = +1.196 C\C F - 0.592 C A Cj + 1.377 C% 

- 1.130 C A C F T F n f - 0.015 C 2 F T F n f - 0.009 CpTjn) , (B.18) 

c 0 s ; 2 (3) = - 2.195 C\C F - 0.649 C A Cj - 1.278 C% 

+ 1.623 C A C F T F nf + 0.244 C 2 F T F n f + 0.025 C F Tjn) , (B.19) 

c 0 s ;i 3) = + 1.058 C\C F + 1.750 C A C 2 F + 0.352 C% 

- 0.635 C A C F T F n f - 0.531 C 2 F T F n f + 0.069 CpTjn) , (B.20) 

with nf = nh + n F 

In addition to the listed coefficient Cj) 3 ^ we require the corresponding coef¬ 
ficient C'^oLg hi the low-energy expansion of the flavour diagonal axial-vector 
correlator in order to compute the correction to the p parameter. Since the 
pole parts of these two coefficients have to cancel, we can decompose the latter 
coefficient as 



4 -^,(3) 

-i,i 


(B.21) 


with the coefficients cjyf as in equations IB.4MB.7l The remaining finite part 


is given by 



+ 2.484 C\C F - 8.319 C A Cj + 16.954 C% 


fin 


- 5.300 C A C F T F n h + 2.759 C 2 F T F n h 

- 1.598 C A C F T F m - 4.210 C 2 F T F m 

- 0.247 C F Tjn 2 h + 1.585 C F Tjn h m + 1.492 CpTjnl , (B.22) 


11 







References 


References 

[1] Y. Schroder, M. Steinhauser, Four-loop singlet contribution to the p 
parameter, Phys. Lett. B622 (2005) 124-130. arXiv:hep-ph/0504055, 
doi:10.1016/j.physletb.2005.06.085. 

[2] K. Chetyrkin, M. Faisst, J. H. Kiihn, P. Maierhofer, C. Sturm, Four-loop 
QCD corrections to the p parameter, Phys. Rev. Lett. 97 (2006) 102003. 
arXiv:hep-ph/0605201, doi:10.1103/PhysRevLett.97.102003. 

[3] R. Boughezal, M. Czakon, Single scale tadpoles and O(GFm^al) 
corrections to the p parameter, Nucl. Phys. B755 (2006) 221-238. 
arXiv:hep-ph/0606232, doi:10.1016/j.nuclphysb.2006.08.007. 

[4] A. Penin, M. Steinhauser, Heavy light meson decay constant from QCD 
sum rules in three loop approximation, Phys.Rev. D65 (2002) 054006. 
arXiv:hep-ph/0108110, doi : 10.1103/PhysRevD. 65.054006. 

[5] M. Jamin, B. O. Lange, fs and fs e from QCD sum rules, 

Phys.Rev. D65 (2002) 056005. arXiv:hep-ph/0108135, 

doi:10.1103/PhysRevD.65.056005. 

[6] I. Allison, et ah, High-precision charm-quark mass from current-current 
correlators in lattice and continuum QCD, Phys. Rev. D78 (2008) 054513. 
arXiv:0805.2999, doi:10.1103/PhysRevD.78.054513. 

[7] C. McNeile, C. Davies, E. Follana, K. Hornbostel, G. Lepage, High- 
precision c and b masses, and QCD coupling from current-current cor¬ 
relators in lattice and continuum QCD, Phys. Rev. D82 (2010) 034512. 
arXiv:1004.4285, doi:10.1103/PhysRevD.82.034512. 

[8] K. G. Chetyrkin, J. H. Kiihn, C. Sturm, Four-loop moments 
of the heavy quark vacuum polarization function in perturbative 
QCD, Eur. Phys. J. C48 (2006) 107-110. arXiv: hep-ph/0604234, 
doi:10.1140/epj c/s2006-02610-y. 

[9] R. Boughezal, M. Czakon, T. Schutzmeier, Charm and bottom 
quark masses from perturbative QCD, Phys. Rev. D74 (2006) 074006. 
arXiv:hep-ph/0605023, doi:10.1103/PhysRevD.74.074006. 

[10] C. Sturm, Moments of heavy quark current correlators at four-loop or¬ 
der in perturbative QCD, JHEP 09 (2008) 075. arXiv:0805.3358, 
doi:10.1088/1126-6708/2008/09/075. 

[11] A. Maier, P. Maierhofer, P. Marquard, The second physical moment of the 
heavy quark vector correlator at 0(a)!), Phys. Lett. B669 (2008) 88-91. 

arXiv:0806.3405, doi:10.1016/j.physletb.2008.09.041. 


12 





[12] A. H. Hoang, V. Mateu, S. Mohammad Zebarjad, Heavy quark vacuum 
polarization function at 0(a 2 ) and 0(a^), Nucl. Phys. B813 (2009) 349- 
369. arXiv:0807 . 4173, doi:10.1016/j.nuclphysb.2008.12.005. 

[13] A. Maier, P. Maierhofer, P. Marquard, A. V. Smirnov, Low energy moments 
of heavy quark current correlators at four loops, Nucl. Phys. B824 (2010) 

1 18. arXiv:0907.2117, doi:10.1016/j.nuclphysb.2009.08.011 

[14] Y. Kiyo, A. Maier, P. Maierhofer, P. Marquard, Reconstruction of heavy 
quark current correlators at O(a^), Nucl. Phys. B823 (2009) 269-287. 
arXiv:0907.2120, doi:10.1016/j.nuclphysb.2009.08.010. 

[15] V. Novikov, L. Okun, M. A. Shifman, A. Vainshtein, M. Voloshin, 

V. I. Zakharov, Sum rules for charmonium and charmed mesons decay 
rates in Quantum Chromodynamics, Phys. Rev. Lett. 38 (1977) 626. 
doi:10.1103/PhysRevLett.38.626 

[16] V. Novikov, L. Okun, M. A. Shifman, A. Vainshtein, M. Voloshin, 

M. Voloshin, Charmonium and gluons: Basic experimental 

facts and theoretical introduction, Phys. Rept. 41 (1978) 1-133. 
doi:10.1016/0370-1573(78)90120-5. 

[17] K. Chetyrkin, J. Kuhn, A. Maier, P. Maierhofer, P. Marquard, 

M. Steinhauser, C. Sturm, Charm and bottom quark masses: 

An update, Phys. Rev. D80 (2009) 074010. arXiv: 0907.2110, 
doi:10.1103/PhysRevD.80.074010. 

[18] J. Koponen, et ah, Heavy-light current-current correlators, PoS LAT- 
TICE2010 (2010) 231. arXiv: 1011.1208. 

[19] A. Djouadi, P. Gambino, Electroweak gauge bosons selfenergies: 

Complete QCD corrections, Phys. Rev. D49 (1994) 3499-3511. 
arXiv:hep-ph/9309298, doi:10.1103/PhysRevD.49.3499. Erratum 
Phys. Rev. D53 (1996) 4111, doi:10.1103/PhysRevD.53.4111 

[20] A. Djouadi, P. Gambino, QCD corrections to higgs boson selfener¬ 
gies and fermionic decay widths, Phys. Rev. D51 (1995) 218-228. 
arXiv:hep-ph/9406431. doi:10.1103/PhysRevD.51.218. Erratum Phys. 
Rev. D53 (1996) 4111, doi : 10.1103/PhysRevD . 53.4111.2.. 

[21] K. Chetyrkin, M. Steinhauser, Three loop nondiagonal current cor¬ 
relators in QCD and NLO corrections to single top quark pro¬ 
duction, Phys. Lett. B502 (2001) 104-114. arXiv :hep-ph/0012002, 
doi:10.1016/S0370-2693(01)00179-4, 

[22] K. Chetyrkin, M. Steinhauser, Heavy - light current correlators at 
order a 2 in QCD and HQET, Eur. Phys. J. C21 (2001) 319-338. 
arXiv:hep-ph/0108017, doi:10.1007/sl00520100744. 


13 




[23] A. Maier, P. Marquard, Low- and high-energy expansion of heavy-quark 
correlators at next-to-next-to-leading order, Nucl. Phys. B859 (2012) 1-12. 
arXiv:1110.5581. doi:10.1016/j.miclphysb.2012.01.021. 

[24] J. Hoff, M. Steinhauser, Moments of heavy-light current correlators up 
to three loops, Nucl. Phys. B849 (2011) 610-627. arXiv: 1103.1481. 
doi:10.1016/j.nuclphysb.2011.04.007 

[25] J. Grigo, J. Hoff, P. Marquard, M. Steinhauser, Moments of heavy 

quark correlators with two masses: exact mass dependence to 

three loops, Nucl. Phys. B864 (2012) 580-596. arXiv: 1206.3418. 
doi:10.1016/j.nuclphysb.2012.07.007 

[26] P. Baikov, K. Chetyrkin, J. H. Kuhn, Vacuum polarization in pQCD: 
First complete 0(a^) result, Nucl. Phys. Proc. Suppl. 135 (2004) 243-246. 

doi:10.1016/j.nuclphysbps.2004.09.013. 

[27] P. Baikov, K. Chetyrkin, J. Kilim, R(s) and hadronic r-decays in order 
O(af): Technical aspects, Nucl. Phys. Proc. Suppl. 189 (2009) 49-53. 

arXiv:0906.2987. doi:10.1016/j.nuclphysbps.2009.03.010 

[28] K. Chetyrkin, A. Maier, Massless correlators of vector, scalar and ten¬ 
sor currents in position space at orders af and of: Explicit ana¬ 
lytical results, Nucl. Phys. B844 (2011) 266-288. arXiv: 1010.1145, 
doi:10.1016/j.nuclphysb.2010.11.007 

[29] P. Nogueira, Automatic Feynman graph generation, J. Comput. Phys. 105 
(1993) 279-289. doi:10.1006/jcph.1993.1074. 

[30] J. Vermaseren, New features of FORM. arXiv :math-ph/0010025. 

[31] M. Tentyukov, J. A. M. Vermaseren, The multithreaded version of FORM, 
Comput. Phys. Commun. 181 (2010) 1419-1427. arXiv: hep-ph/0702279, 
doi:10.1016/j.cpc.2010.04.009. 

[32] T. van Ritbergen, A. Schellekens, J. Vermaseren, Group theory fac¬ 
tors for Feynman diagrams, Int. J. Mod. Phys. A14 (1999) 41-96. 
arXiv:hep-ph/9802376, doi:10.1142/S0217751X99000038. 

[33] C. W. Bauer, A. Frink, R. Kreckel, Introduction to the GiNaC framework 
for symbolic computation within the C++ programming language, J. Symb. 
Comput. 33 (2000) 1. arXiv: cs/0004015 

[34] R. Lewis., Fermat 

URL https://home.bway.net/lewis/ 

[35] P. Marquard, D. Seidel, Crusher. 

[36] S. Laporta, High precision calculation of multiloop feynman integrals 
by difference equations, Int. J. Mod. Phys. A15 (2000) 5087-5159. 
arXiv:hep-ph/0102033, doi:10.1016/S0217-751X(00)00215-7. 


14 




[37] Y. Schroder, A. Vuorinen, High-precision epsilon expansions of 
single-mass-scale four-loop vacuum bubbles, JHEP 0506 (2005) 051. 
arXiv:hep-ph/0503209, doi:10.1088/1126-6708/2005/06/051. 

[38] M. Faisst, P. Maierhofer, C. Sturm, Standard and e-finite mas¬ 
ter integrals for the p-parameter, Nucl. Plrys. B766 (2007) 246-268. 

arXiv:hep-ph/0611244, doi:10.1016/j.nuclph.ysb.2006.12.014. 

[39] P. Maierhofer, Four-loop vacuum diagrams [in German], Diploma Thesis 
(Karlsruhe 2006). 

[40] M. Veltman, Limit on mass differences in the Weinberg Model, Nucl. Phys. 
B123 (1977) 89. doi:10.1016/0550-3213(77)90342-X. 

[41] A. Pak, The toolbox of modern multi-loop calculations: novel analytic 
and semi-analytic techniques, J. Phys. Conf. Ser. 368 (2012) 012049. 
arXiv:1111.0868. doi:10.1088/1742-6596/368/1/012049. 


15 


